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AN INITIAL INVESTIGATION INTO METHODS OF COMPUTING 
TRANSONIC AERODYNAMIC SENSITIVITY COEFFICIENTS 


I. Introduction 

This report covers the period from July 1991 thru December 1992. During this 
reporting period, a method based upon the quasi-analytical approach has been developed 
for computing the aerodynamic sensitivity coefficients for three-dimensional wings in 
transonic and subsonic flow. In addition, the method computes for comparision purposes 
the aerodynamic sensitivity coefficients using the finite-difference approach. The 
accuracy and validity of the methods are currently under investigation. 


H. Personnel 

The individuals associated with the project during this reporting period have been 
Dr. Leland A. Carlson, Principal Investigator, and Hesham Elbanna, Graduate Research 
Assistant. Mr. Elbanna has been supported by the project during this period and will use 
the results of this research effort for his Ph.D. dissertation. 


m. Research Progress 

During the past six months, significant progress has been achieved in developing the 
quasi-analytical approach to obtaining aerodynamic sensitvity coefficients about wings in 
transonic flow. As mentioned above, a method and computer program has been 
developed which can compute such coefficients for wings in subsonic and transonic flow. 
In addition, the method also automatically computes for comparison purposes the 
aerodynamic sensitvity coefficients using the finite-difference approach. 

The method consists of several fundamental components. The first is the transonic 
flow solver which is a three-dimensional full potential method using the zebra algorithm 
as developed by Carlson, Weed, and Anderson. The flow solver uses are Cartesian like 
grid and places the surface boundary conditions in the x-y plane. It is a very robust and 
efficient flow solver and should be adequate for the present studies. 

The second portion determines the "analytical" sensitivity derivatives using the 
quasi-analytical approach. The code for this task was developed using Macsyma to 
determine appropriate derivatives, i.e. dR/30's and 5R/3X D 's, where R is the residual, </) 
is the potential, and X D 's are the design variables. The resultant set of algebraic 
equations for the medium grid being considered is 17500 x 17500. These algebraic 
equations are solved using an IBM3090 conjugate gradient solver. For the sensitivity 
coefficient solver, twelve basic design variables and six derived design variables are 
considered. The basic design variables are — ffeestream Mach number, angle of attack, 
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ma ximu m airfoil thickness, maximum airfoil camber, location of maximum camber, 
twist at four locations, wing tip leading edge streamwise coordinate, wing tip trailing 
edge streamwise coordinate, and wing tip spanwise coordinate. From these, the derived 
design variables considered are semi-span, wing area, aspect ratio, taper ratio, leading 
edge sweep angle, and trailing edge sweep angle. 

In addition, to the quasi-analytical approach, the current computer code and method 
also computes sensitivity derivatives using the finite-difference approach. In this section, 
the input design variables are changed automatically and the transonic flow solver is re- 
run to obtain a new solution. The derivatives are then determined using j for each 
grid point in the flow field. From these values, the aerodynamic sensitivity coefficients 
can then be determined. 

The method also contains a section which creates and plots extensive graphical 
output. For example the bC p (x)/dX D distribution and the 3Q/dX D is computed at twenty 
spanwise stations and plotted at user selected stations. Further, spanwise variations of 
section senstivity derivatives are determined and plotted along with predicted C p (x) 
distributions determined using the computed sensitvitiy derivatives. 

It should be noted that currently there are two versions of the quasi-analytical 
sensitivity derivative code. The first computes “bR/d^'s ignoring any dependence of the 
potential flow upwind switch funtion on the (p's. This version is called the Nu = C 
version. The second scheme computes the oR/oO's including the dependence of the 
upwind switching function on the $ values; and it is termed the Nu=f($) version. 

At this point it is recognized that neither approach is perfect or validated. Further, 
to date all calculations have been executed using only single precision arithmetic. Based 
upon previous two-dimensional studies, it is suspected that the finite difference approach 
may require double precision execution in order to yield correct values. Thus, in the 
comparisions which follow, the fact that the quasi-analytical approach and the finite 
difference method yield different values for the sensitivity derivatives does not 
necessarily imply that the quasi-analytical method is in error. 

In the following section, the viewgraphs used in a presentation to NASA Langley 
are reproduced. These chart effectively summarize the current state of the research and 
are indicative of the type of results which can be obtained from the present approach. 
While not definitely established, it is believed that the present Nu=f(0) version has the 
correct behavior and is the better of the two quasi-analytical versions. However, it is 
believed that it still contains some "errors", and this possibility is currently under 
investigation. 
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IV. Future Efforts 


During the next reporting period, work will continue on developing the quasi- 
analytical approach and verifying its usefulness from a proof-of-concept viewpoint. In 
addition, it is planned to prepare and present a paper at the 1992 AIAA Applied 
Aerodynamics Meeting and at the 3rd Pan American Congress of Applied Mechanics. 

V. Technical Monitor 

The technical monitor for this project is Dr. Woodrow Whitlow, Jr., Unsteady 
Aerodynamics Branch, MS 173, NASA Langley Research Center. Dr. Whitlow replaces 
Dr. E. Carson Yates, Jr., who retired recently from the Interdisciplinary Research Office 
at NASA Langley. 
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CHARTS USED IN THE PROGESS REPORT 
PRESENTATION TO 
NASA LANGLEY 
ON 

FEBRUARY 24, 1992 
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Important Points 


This is a report on research in progress. 

The computer codes are still under 
development. 

The answers are not perfect. They may even 
be wrong. 

There is still a lot of work to be 
accomplished. 

Primary objective is investigate the quasi- 
analytical approach to aerodynamic 
sensitivity derivatives, determine methods 
for finding the derivatives using the QA 
approach, and establish its "validity" and 
range of applicability. I.E. PROOF OF 
CONCEPT. 



Accomplishments 


A method has been developed for three- 
dimensional flow to compute aerodynamic 
senstivity coefficients using the quasi- 
analytical approach. 


The method also computes for comparison 
aerodynamic sensitivity coefficients using 
the finite-difference approach. 



Components of Method 


Flow Solver 

3-D full potential method using the 
Zebra algorithm as developed by Carlson, 
Weed, and Anderson. Zebra uses are 
Cartesian like grid and places surface 
boundary conditions in the x-y plane. 

Good Idea ?? Yes and No. 


"Analytical" Sensitivity Derivatives 

Based on quasi-analytical approach. 
Code developed using Macsyma to 
determine appropriate derivatives, 5R/60 's 
and 5R/(3 Xd 's. 

Resultant set of algebraic equations is 
about 17500 x 17500. Solved for 50/5Xd 
values using IBM3090 conjugate gradient 
solver. 



Components of Method 

(continued) 

Finite-Difference Sensitivity Derivatives 

Input design variables changed 
automatically and flow solver run to obtain 
new solution. Derivatives determined using 
^.0 /^Xd for each grid point. 


Graphical Output 

5Cp(x)/6Xd distribution and dCl/dXo 
computed at 20 spanwise stations and 
plotted at selected stations. 

Spanwise variations of section sensitivity 
derivatives determined and plotted. 

Predicted Cp(x) distributions determined 
using sensitivity derivatives and plotted. 

Etc. Etc. Etc. 



Twelve Basic Design Variables 

Freestream Mach Number 
Angle of Attack 
Maximum Airfoil Thickness 
Maximum Camber 
Location of Maximum Camber 
Twist at Four Locations 
Wing Tip Leading Edge Coordinate 
Wing Tip Trailing Edge Coordinate 
Wing Tip Coordinate 


Six Derived Design Variables- 

Semi-span 
Wing area 
Aspect Ratio 
Taper Ratio 
Leading Edge Sweep 

Trailing Edge Sweep 



Comments 


Currently there are two versions of the 
quasi-analytical sensitvity derivative code. 
The first computes 3R/<30's ignoring any 
dependence of the potential flow upwind 
switch function on (£>'s. This is called the Nu 
= C version. The second computes dR/00's 
including the dependence of the upwind 
switching function on 0's. It is termed the 
Nu=f(0) version. 

Neither method is "perfect" or "validated". 

Not all known "changes" have been 
included. 

In the results, you will see solutions from 
various versions of each approach. 



Current Situation 


We are suffering from a bad case 

of 
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Figure 1 — Conditions for Subcritical Test Case 
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Fig. 4B -- Quasi-Analytical and Finite Diff. Basic Sensitivity Derivatives 
M 00 = 0.8, AOA = lo, Nu = f(0) in QA, AX D = 10-3 
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Fig. 5 — Derived Sensitivity Derivatives 
M 00 = 0.8, AOA = 1°, Nu = {((])) in QA, AX D = 10’ 3 
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Fig. 6 -- Conditions for First Supercritical Test Case 
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Fig. 8A — Quasi-Analytical and Finite Diff. Basic Sensitivity Derivatives 
M 00 = 0.8, AOA = 3°, Nu = f (<£)) in QA, AXq = 10' 3 
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Fig. 9 -- Derived Sensitivity Derivatives 
0.8, AOA = 3°, Nu = f((f>) in QA, A X D = 10 ' 3 
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Fig. 10 - Conditions for Second Supercritical Test Case 









Fig. 1 1 -- Cp Distribution for Supercritical Test Case, Cp* = -0. 379 


















aCp/8L.ACpML aCp/aT.ACp/AT aCp/8M.ACp/AM 



20.0 


10.0 


0.0 


10.0 


20.0 





Fig. 12A — Quasi-Analytical and Finite Diff. Basic Sensitivity Derivatives 
M oc = 0.82, AOA = 3°, Nu = f(0) in QA, Md = 10' 3 
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Fig. 12B - Quasi-Analytical and Finite Diff. Basic Sensitivity Derivatives 
M 00 = 0.82, AOA = 3°, Nu = f(0) in QA, 2 $X D = lO 3 
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Fig. 13 -- Derived Sensitivity Derivatives 
M 00 = 0.82, AOA = 3°, Nu = f(<fo in QA, AX D = 1°' 3 
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Fig. 14 -- Conditions for Third Supercritical Test Case 
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Fig. 17 - Spanwise Variation of Sensitivity Derivatives, M 00 - 0.84, 
AHA = in Nu = f((/)) in QA, AX D =10'3 
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Fig. 18 -- Derived Sensitivity Derivatives 
M 00 = 0.84, AOA = 3°, Nu = f(0) in QA, &X D = 10-3 
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Fig. 20A -- Quasi-Analytical and Finite Diff. Basic Sensitivity Derivatives 
M no = 0.84, AOA = 3°, Nu = C in QA, AX D = 10 * 
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Fig. 22 -- Derived Sensitivity Derivatives 
= 0.84, AOA = 3°, Nu = C in QA, &X D = 10'3 
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Fig. 23 -- Spanwise Variation of Derived Sensitivity Derivatives, 
M of) = 0.84, AOA = 3°, Nu = C in QA, AX D = 10* 3 













NOTE 


Comparison of Execution Times 

Derivatives by Finite Differences = 
9.8 

Derivatives by QA Nu = C version 
= 3.9 

Derivatives by QA Nu = f(<j£)) 
version = 5.8 



Question 

Which version is "correct"? 
Nu = f(0) 
or 


Nu = C 


To determine, consider some of our 
previous 2-D results. 
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Figs. 25 -- Two-Dimensional Results at M 00 = 0.82 
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^ Figs. 25 -- Two-Dimensional Results at M 00 - 0.82 
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Figs. 25 — Two-Dimensional Results at M 00 = 0.82 
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Figs. 26 -- Two-Dimensional Results at M 00 = 0.84 
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Figs. 26 -- Two-Dimensional Results at M 00 0.84 
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Figs. 26 — Two-Dimensional Results at M 00 = 0.84 
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Conclusions 


The Nu =f(< £)) version of the QA 
method has the same behavior and 
trends as our previous 2-D results. 

Differences from QA results and 
FD results are larger in present 3-D 
cases than in previous 2-D cases. 
While, current results are all single 
precision, these differences 
indicate that "errors" probably still 
exist in the method. 

Question 

What is the influence of the size of 
AXd on the finite difference 
results? 
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Fig. 27A - Effect of &M 0O Size on Fin. Diff. Denvatives 
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Fig. 27B - Effect of aM 00 Size on Fin. Diff. Derivatives 







th = 0.84 
<V = 3° 

* 1 -=. 0 , 564 - 



r;« __ PTIVv't nf A M Size nn ftCi/flM hv Fill. Diff. 


3Cp/3ot,ACp/Aa 3Cp/3a,ACp/Aa 











dCp/3oc,ACp/A a 


7](12)= 0.564 

- 10 ' L 










5.00 


4.50 


4.00 


5.50 


5.00 


2.50 


2.00 

0.00001 



IV/dOV'ie/dDP IV/dDV'ie/dDe 







Fig. 3 1 A -- Effect ofAT Size on Fin. Diff. Derivatives 
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Fig. 3 IB — Effect of AT Size on Fin. Diff. Derivatives 
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Fig. 33A - Effect of AC Size on Fin. Diff. Derivatives 
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Pig. 33B - Effect of LC Size on Fin. Diff. Derivatives 
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Fig. 35 A - Effect of AL Size on Fin. Diff. Derivatives 
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Question 


What is the ability of the present 
aerodynamic sensitvity coefficients to 
predict C p (x) distributions away from the 
nominal point? i.e. 

C p (x,y,Xo+ilXD) = C p (x,y,Xo) + 
(dC p /dX D )x,y * £XD 

To be fair, the sensitivity derivative in the 
finite difference case must be evaluated at a 
(hX d)i different from aXd- Also, the aXd's 
must be reasonable. For the following, in 
the FD cases, the derivatives were computed 
using 0.001 increments. 




Fig. 37 - Comparison of Cp Predictions using QA Derivatives at Nominal 
Point with Actual Results, Nu = f(Q) in QA. 
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Fig. 39B — Comparison of Cp Predictions Using QA with Nu variable, QA 
with Nu constant, and Fin. Diff. Derivatives at Nominal Point with Actual 

Results. 









Conclusions 


Believe Nu = f((j)) version is better of the two. However, it 
still needs improvement. Plan to have a reasonable version 
by summer. 

Currently having Macsyma computer problems. 

If know AXd needed, a finite difference approach will 
obviously yield correct results. If don't know Z\Xd, QA 
approach will give reasonable estimates of aerodynamic 
sensitivity derivatives. 

Present code needs some further generalization. Depends 
upon funding situation, new student learning method, etc. 

Might be desirable to "repeat" process with extended small 
disturbance equations. Matrix elements would be easier to 
determine. Result might be more general. 

Will present a paper at 1992 Applied Aerodynamics 
Meeting. Hope to also present one in Sao Paulo in January 
1993. 



